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1. Introduction 



Photoproduction of heavy flavours has been studied extensively in fixed target 
experiments (ref. [|T|). The measured total cross sections for charm photoproduction 
turn out to be in reasonable agreement with theoretical expectations (see ref. 0). 
At HERA new opportunities are available. On one hand, one can extend in en- 
ergy previous charm production studies, from center-of-mass energies around 30 GeV 
to energies up to 300 GeV. Some prehminar results have already been presented 
in ref. [^. Small-x effects in these high energy regimes may not be negligible for 
charm production, and it will be interesting to see whether they are important at 
HERA. Bottom production will also be observable at HERA. The larger value of 
the bottom quark mass makes QCD predictions more reliable. Furthermore, bottom 
hadropro duct ion has been extensively studied at hadron colliders, and the compari- 
son with QCD predictions, although qualitatively good, presents some problems (see 
refs. Theoretical arguments indicate that perturbation theory should de- 

scribe photoproduction of heavy flavours more reliably than hadropro duct ion. This 
is also supported by experimental evidence in fixed-target charm production. HERA 
has therefore the opportunity to clarify some open questions in bottom production. 

In this paper, we discuss the total cross sections for charm and bottom photo- 
production, and for electroproduction in the Weizsacker- Williams approximation, in 
the HERA energy regime. Our analysis is based on the next-to-leading order calcula- 
tion of heavy-quark photoproduction and hadroproduction cross sections performed in 
refs. . Some results in heavy-quark electroproduction at HERA have already been 
presented in ref. [|^. The computer program we used for this work was developed 



in refs. |Tl|,|T2|, and is available upon request. We will not consider electroproduction 
cross sections for large photon virtuality, which has been recently discussed in ref. 0. 

Our work is organized as follows. In section 2 we discuss charm photoproduc- 
tion, focussing especially on the theoretical uncertainties which are specific of HERA 
physics. In section 3 we give our results for bottom photoproduction. In section 4 we 
present our results for charm and bottom production in ep collisions, and in section 5 
we give our conclusions. 



-2- 



2. Charm production 

In charm photoproduction at HERA, because of the large center-of-mass energies 
available, the new problem arises (which was not present in fixed target configura- 
tions) of the presence of a large radiative effect, proportional to aemoi log(5'/4m^), 
where S is the squared center-of-mass energy of the photon-hadron system, and rric 
is the charm quark mass. This effect is the first manifestation of a whole tower of 
corrections behaving like ae^ag{aslog{S/4m^))^, arising from the small-x region. In 
the high energy limit, log(S'/4m^) may become of order 1, and one should then 
worry about resumming all these contributions. A formalism for the resummation has 
been studied in refs. ||T^. In the present work, we will try to assess the importance 
of these effects at HERA. Following the work of ref. in fig. we show various 
contributions to the total pointlike photoproduction cross section, histogrammed as 
a function of the variable r, defined as the ratio of the partonic versus the hadronic 
center-of-mass squared energies s/S. 




Figure 1: The ©(aemOs) 79, C'(aemas) 19 '^"■^ 0(acma|) 79 contributions 
to the total charm photoproduction cross section, histogrammed as a function 
of T = s/S. The proton factorization and renormalization scales have been 
taken equal to ruc in the short- distance cross sections. 



In the case of the pointhke component of photoproduction, this ratio is simply 
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equal to x, the fraction of the hadron momentum carried by the parton. We see 
that in the next-to-leading 'yg contribution there is a sizeable hump away from the 
threshold region. The same behaviour is observed in the 'jq contribution. The hump 
arises from the fact that the next-to-leading parton cross section goes to a constant 
for large s (see figs. 2 and 3 of ref. Therefore, for a hypothetical gluon density of 
the form G{x) = 1/x, the next-to- leading 'yg contribution to the total cross section 
would be given by 

o-W ^(a;) dx = r a^dlogx. (2.1) 

Observe that, were it not for the fact that the gluon density has some extra damping 
factor as x 1, the dashed histogram in fig. ^ would go to a constant as a; 1. In 
the very high energy limit, the hump will therefore become a plateau, and its width 
would be roughly equal to log(S'/4m^). One can show that this problem does not stop 
at this order: at order aem^s the partonic cross section grows like (log s/4m^)"~^, 
and this would give a contribution to the total cross section raising like (log(5'/4m^))'^ 
with the energy. Notice that fig. |^ was obtained using the MRSA parton densities 



T4| , which have a relatively singular gluon distribution at small x. This type of 
parton density emphasizes the threshold region. We therefore expect that when using 
parton densities with a more regular behaviour at small x, the hump should be more 
pronounced. This is in fact the shown in fig. ^, obtained using the parton 

densities HMRS BQ |15[ 

The size of the hump relative to the total cross section can be taken as a measure 
of the importance of the small-x regime. In fig. |^ we plot the contribution to the 
total cross section coming from the hump, relative to the total and to the (9(agaem) 
contribution. The hump contribution is defined as the integral of the (9(a|aem) 
term in da/dr, from its minimum point up to r = 1. We are thus led to conclude 
that resummation effects could be important at the typical center-of-mass energies 
explored by HERA. Furthermore, flatter parton densities lead to larger small-x effects. 
A similar pattern was found in the context of hadroproduction of bottom quarks, in 
ref. , where the all-order resummation of small- a; effects was actually performed for 
h production at the Tevatron. Lacking a calculation of the effects of resummation in 



^In this case the parton density is frozen to its value at the minimum ahowed by the 
parametrization. 
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Figure 2: As in fig. ||, using HMRS B parton densities. 
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Figure 3: The contribution of the hump to the total cross section relative 
to the 0(a|aem) contribution (upper) and to the total (lower) for the MRSA 
and HMRS B proton parton densities. 
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photoproduction, we will simply assume that it may give an enhancement equal to the 
square of the size of the hump relative to the total cross section, a pattern suggested 
by the results of ref. |16]. From fig. ^ we can then conclude that our calculation 



may underestimate the total charm cross section by 20% to 40%, depending upon the 
parton density choice. 

Sizeable small-x effects are also present in the hadronic component of the charm 
cross section. In this case, we find that the relative importance of the hump contri- 
bution varies between 20% and 45% when the center-of-mass energy changes from 
50 to 300 GeV (this result was obtained using the MRSA and GRV-HO ^ parton 
densities in the proton and the photon respectively). The corresponding uncertainty 
on the final result is considerably smaller than other uncertainties. 

Other important inputs in the computation of the total cross section are the value 
of Af^ and the gluon density in the proton. The world average for as{M^) taken 



from the Particle Data Book, ref. |18], is as{M^) = 0.117 ± 0.005, corresponding 



to Af^ = 215 + 65 — 55MeV (the value of A^^ quoted in ref. [|T8| differs from ours 
because it actually refers to the 3- loop formula for as, while we are quoting the 2- 
loop value). On the other hand, the preferred value of Af^ used in the latest parton 
density parametrizations of the CTEQ group is 158 MeV (see ref. [|19|), while the 
MRSA parametrization uses 151 MeV (see ref. Jl^)- These values are below the 
lowest extreme of the world average range. This is due to the well known fact that 
the determinations of ag from e'^e~ data tend to be higher than those from deep- 
inelastic scattering. In order to assess the error in our predictions coming from the 
uncertainties on Af^, we will therefore consider also the parton density set CTEQ2ML, 
which was fitted with A^^ = 220 MeV, a value close to the world average. In summary, 
our range for Af^ is 151 MeV< Af^ < 220 MeV, which does not cover the full 
uncertainty on the present knowledge of a^. We are forced to do so, since larger 
values of Af^ are too far from deep inelastic scattering determinations, and it would 
be difficult to assess their correlation with the gluon density in the proton. 

Several uncertainties come into play in the extraction of the gluon density. Its 
normalization is strongly constrained by the momentum sum rule, while its shape, 
at relatively moderate value of x, is constrained by direct photon production data. 
Deep inelastic scattering data constrain the gluon density for relatively small values 
of X, where it affects the evolution of the singlet parton density. At very small 
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values of x, tlie gluon density is in fact not yet known with precision. Both the 
CTEQ and MRS groups assume that its small-x behaviour is the same as that of 
F2, which is measured at HERA down to values of x of the order of 10~^. On the 
other hand, a measurement of the gluon density at Q"^ = 7 GeV^ has recently been 
performed by the ZEUS collaborationil. We report in fig. ^ the gluon densities 
given by the sets MRSA, MRSD-' and CTEQ2MF, together with the two extreme 
small-x parametrizations obtained by the ZEUS collaboration. The sets MRSD-' and 

40 
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X 

Figure 4: The gluon densities in MRSA, MRSD-' and CTEQ2MF, together 
with the upper and lower parametrization of the ZEUS measurement. 

CTEQ2MF differ essentially in the assumed small-x behaviour. They are consistent 
with the two extreme parametrizations of the ZEUS data, and we will therefore use 
them as the two extreme possibilities. 

Our results for the pointlike component of the charm cross section are shown 
in table |1| and in figs. ^ and ^. The strong sensitivity of the cross section to the 
charm quark mass is quite apparent. Scale sensitivity is also important, giving an 
uncertainty of a factor of 2 for the low mass value. From fig. |^ we also see that by 
comparing the cross section measured at HERA with lower energy measurements it 
should be possible to distinguish between different small-x behaviour of the structure 
functions. 

The hadronic component of the charm cross section is sensitive to the behaviour 
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Proton PDF set MRSA 


30 GeV 


2.582 


2.104 


1.712 


1.321 


1.090 


0.907 


0.735 


0.611 


0.515 


100 GeV 


6.558 


5.056 


4.057 


3.415 


2.826 


2.356 


2.009 


1.721 


1.467 


280 GeV 


13.89 


10.13 


8.008 


7.314 


5.858 


4.831 


4.389 


3.685 


3.120 


Proton PDF set CTEQ2MF 


30 GeV 


2.651 


2.264 


1.885 


1.455 


1.216 


1.023 


0.827 


0.688 


0.583 


100 GcV 


4.371 


3.856 


3.236 


2.699 


2.425 


2.088 


1.761 


1.593 


1.393 


280 GeV 


7.292 


5.681 


4.607 


4.504 


3.827 


3.237 


3.035 


2.676 


2.317 


Proton PDF set MRSD-' 


30 GcV 


3.136 


2.363 


1.885 


1.422 


1.136 


0.936 


0.748 


0.612 


0.513 


100 GeV 


11.75 


8.149 


6.347 


5.106 


3.951 


3.220 


2.690 


2.199 


1.843 


280 GeV 


34.16 


23.06 


17.81 


14.79 


11.17 


9.033 


7.794 


6.247 


5.196 


Proton PDF set CTEQ2M 


30 GeV 


2.696 


2.182 


1.785 


1.341 


1.107 


0.926 


0.734 


0.612 


0.519 


100 GeV 


6.618 


5.240 


4.262 


3.473 


2.917 


2.454 


2.051 


1.773 


1.522 


280 GeV 


13.46 


10.17 


8.157 


7.207 


5.898 


4.919 


4.374 


3.728 


3.184 


Proton PDF set CTEQ2ML 


30 GeV 


3.500 


2.435 


1.829 


1.705 


1.242 


0.974 


0.910 


0.688 


0.553 


100 GeV 


8.818 


5.625 


4.171 


4.270 


3.169 


2.497 


2.451 


1.942 


1.583 


280 GeV 


18.90 


10.52 


7.618 


8.844 


6.205 


4.823 


5.145 


3.972 


3.213 



Table 1: Pointlike component of the total charm cross sections (fib) in jp 
collisions, for various input parameters, as indicated. 
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Figure 5: Pointlike component of the total charm cross section versus the 
center- of -mass energy in 7p collisions. The bands are obtained by varying the 
renormalization scale between mc/2 and 2mc- The proton parton density set 
is MRS A. 
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Figure 6: Pointlike component of the total charm cross section versus the 
center- of -mass energy in 'yp collisions for the three sets of parton densities 
MRSA, MRSD- and CTEQ2MF. The bands are obtained by varying the 
renormalization scale between micjl and 2mc- 
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of the gluon density in the photon at relatively small values of x, a region not yet 
well explored experimentally. In fig. |^ we show the hadronic component cross section 
histogrammed as a function of (the partonic x in the photon). We see that the 




Figure 7: The distribution of the hadronic component for charm produc- 
tion in jp collisions. 



dominant values of x^ are between 10 ^ and 10 ^. Results from jet photoproduction 
at HERA0 give some indication on the gluon density in the photon, whose extraction 
depends upon the value of the quark densities in the photon in a region where they are 
not directly probed. These data, in spite of the uncertainties, seem to indicate that the 
gluon density is between the one in the GRV-HO set and the LACl set p2[ . We 
note that the available experimental information on the gluon density in the photon 
do not allow to distinguish among different next-to-leading order parametrizations 



23| of the parton distributions. Furthermore, the gluon density in the more recent 
leading order parametrization of ref. is within the range spanned by the GRV-HO 
and the LACl sets. We will therefore take the sets GRV-HO and LACl as the two 
extremes. For the proton densities we will use the MRS A set. The relevant results are 
reported in table ^ and in fig. ^. In fig. ^ we plot the cross section without hadronic 
component, and with the hadronic component evaluated with the GRV-HO and LACl 
sets of photon parton densities. We conclude that with the present uncertainty it is 
quite possible that the hadronic component of the photon dominates the cross section 
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at typical HERA energies. The pointlike and hadronic components of the cross section 
may be separated experimentally, at least in a first approximation, since in the case 
of the hadronic component a sizeable fraction of the photon momentum is lost into 
hadronic fragments. Therefore, at HERA it may be possible to distinguish between 
the two sets of parton densities in the photon we have considered. 
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Figure 8: Total charm cross section versus the center- of-mass energy in 
collisions, with and without the inclusion of the hadronic component, computed 
with the two sets of photon parton densities LACl and GRV-HO. The hands 
are obtained by varying the renormalization scale between mc/2 and 2mc- 
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3. Bottom production 

In the case of bottom production, all the uncertainties we have discussed for 
charm are strongly reduced. The small-x problem is much less dramatic. We have 
performed a study of the bottom cross section versus r, analogous to the one presented 
for charm in the previous section, and found that the size of the hump contribution 
is of the order of 20% ot the total cross section. The sensitivity of the pointlike cross 
section to the other parameters is given in table ^. Our central values of the cross 
section, in the column marked with DEF, are obtained with the renormalization 
and factorization scales set equal to mb, and mfe = 4.75 GeV. If we vary all the 
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rric = 1.2 GeV 


= 1.5 GeV 


= 1.8 GeV 


EcM f^R 


mc/2 


rric 


2mc 


mc/2 


rric 


2mc 


mc/2 


rric 


2mc 


Photon PDF set GRV-HO 


30 GeV 


0.793 


0.297 


0.165 


0.206 


0.099 


0.060 


0.072 


0.040 


0.026 


100 GeV 


4.572 


1.680 


0.924 


1.438 


0.673 


0.405 


0.589 


0.314 


0.200 


280 GeV 


15.42 


5.492 


2.989 


5.296 


2.435 


1.454 


2.348 


1.239 


0.784 


Photon PDF set LACl 


30 GeV 


3.350 


0.995 


0.503 


0.540 


0.213 


0.119 


0.125 


0.059 


0.036 


100 GeV 


35.05 


11.69 


6.207 


8.975 


3.794 


2.187 


3.052 


1.469 


0.893 


280 GeV 


140.3 


48.77 


26.30 


42.45 


18.84 


11.09 


16.78 


8.494 


5.283 



Table 2: Hadronic component of the total charm cross sections (fib) in 
7P collisions, for various input parameters, as indicated. The proton parton 
density set is MRS A. 



parameters together in the direction that makes the cross section larger or smaller we 
get a spread for the bottom cross section which goes from 16 nh to 35 nh at a center 
of mass energy of 100 GeV, and from 41 nh to 101 nh at 280 GeV. The results for 
the hadronic component are collected in table ^. The default values correspond to 
fin = fip = fi^ = mb = 4.75 GeV. We see that in this case the effect of the hadronic 
component is less dramatic. According to the LACl set, the hadronic contribution 
could be a substantial fraction of the total, while the GRV-HO set gives a much 
smaller result. 

4. Electroproduction of heavy quarks 

We present here our predictions for total charm and bottom cross sections in 
ep collisions. The total electroproduction cross section in the Weizsacker- Williams 
approximation is given by 

(Tep{S) = / dx f!f\x,fi^^n^)a^p{xS). (4.1) 
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nib (GeV) 
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mh/2 




7715/2 




4.5 
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Proton PDF set MRSA 


100 GeV 


22.90 


26.00 


20.36 


21.87 


23.63 


27.41 


19.28 


280 GeV 


64.30 


75.44 


56.09 


57.81 


66.95 


75.31 


55.32 


Proton PDF set CTEQ2MF 


100 GeV 


24.91 


27.93 


22.33 


23.80 


25.73 


29.70 


21.02 


280 GeV 


62.91 


73.08 


55.21 


56.27 


65.80 


72.76 


51.71 


Proton PDF set MRSD-' 


100 GeV 


23.03 


26.19 


20.15 


22.18 


23.63 


27.75 


19.26 


280 GeV 


74.76 


87.24 


65.43 


69.22 


76.93 


88.76 


63.52 


Proton PDF set GTEQ2M 


100 GeV 


22.96 


25.92 


20.50 


22.04 


23.54 


27.50 


19.31 


280 GeV 


65.33 


75.96 


57.32 


59.43 


67.51 


76.40 


56.23 


Proton PDF set CTEQ2ML 


100 GeV 


25.07 


29.53 


21.76 


23.38 


25.92 


30.03 


21.07 


280 GeV 


71.04 


86.84 


60.42 


61.74 


74.39 


82.99 


61.20 



Table 3: Pointlike component of the total bottom cross sections (nh) in ^yp 
collisions, for various input parameters, as indicated. 
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nib (GeV) 


EcM 


DEF 


mb/2 


2mb 


mb/2 


2mfe 


mb/2 


2mb 


4.5 


5 


Photon PDF set GRV-HO 


100 GeV 


4.244 


5.589 


3.290 


4.246 


4.258 


3.883 


4.525 


5.444 


3.347 


280 GeV 


27.63 


38.67 


20.74 


25.68 


29.13 


25.39 


28.68 


34.29 


22.48 


Photon PDF set LACl 


100 GeV 


6.693 


9.264 


5.053 


7.157 


6.394 


7.161 


6.348 


9.309 


4.890 


280 GeV 


98.45 


137.2 


74.05 


96.07 


100.6 


106.4 


92.18 


128.7 


76.15 



Table 4: Hadronic component of the total bottom cross sections (nh) in 
7P collisions, for various input parameters, as indicated. The proton parton 
density set is MRS A. 



In the following, unless otherwise specified, we will refer to the case when no anti-tag 
condition is imposed on the outgoing electron. The function /^^^ is given in this case 
by 

/f = ^l±ll^iog^fl24^. (4.2) 



The scale fiww, according to ref. ||2^, will be taken equal to the quark mass for the 
pointlike component. 

The choice of fi^vv in the case of the hadronic component is a debatable issue. 
Intuitive reasoning leads to conclude that this scale should not be larger than a typical 
hadronic scale, because a strongly off-shell photon would have smaller probability of 



mixing with hadronic resonances. Some authors (see refs. p6| , p7| , p8| ) suggest instead 
that this scale should be related to the hardness of the process, and should therefore 
be larger for bottom production than for charm production. In practice, in the charm 
case, changing fi^vw from 1 GeV up to the quark mass has only a few percent effect on 
the size of the hadronic component for charm production, while for bottom production 
the effect could be as large as 20%. In both cases, as we will see, these effects are, 
at present, much smaller than other sources of uncertainty. In the following we will 
choose fiww = 1 GeV for charm and ^ww = 3 GeV for bottom. We also notice that, 
in the case when an anti-tagging condition is imposed, this problem is not relevant 
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mc = 1.2 GeV 


nic = 1.5 GeV 


rric = 1.8 GeV 


fJ'R 


mj2 


rric 


2mc 


mc/2 


rric 


2mc 


mc/2 


rric 


2mc 


MRSA 


1.478 


1.113 


0.888 


0.749 


0.605 


0.500 


0.431 


0.361 


0.305 


CTEQ2MF 


1.130 


0.923 


0.761 


0.630 


0.535 


0.453 


0.385 


0.334 


0.288 


MRSD-' 


2.472 


1.713 


1.333 


1.087 


0.833 


0.677 


0.574 


0.463 


0.387 


CTEQ2ML 


2.078 


1.250 


0.918 


0.956 


0.678 


0.529 


0.530 


0.406 


0.327 



Table 5: Pointlike component of the total charm cross sections (fj,b) in ep 
collisions, at Eqm = 314 GeV, for various input parameters, as indicated. 



at all. 

We present in tables ^ and ^ our predictions for the total charm cross sections, 
for the pointlike and the hadronic component respectively, at -EcAf = 314 GeV. The 
proton and photon parton densities are chosen as in the case of photoproduction. We 
do not consider here the CTEQ2M set, which gives essentially the same results as our 
default set MRSA. We observe that the most important source of uncertainty is the 
charm mass. The ratio bewteen the extreme choices rric = 1-2 and 1.8 GeV can be 
as large as 4.3 for the pointlike component and 10 for the hadronic component, and 
becomes larger for smaller values of the renormalization scale. The renormalization 
scale variation has a strong effect on the hadronic component (a factor up to 5.5 for 
the lowest value of the charm mass), and quite a sizeable effect on the pointlike com- 
ponent. As in the photoproduction case, the differences in the hadronic component 
evaluated with GRV-HO and LACl densities are striking. For our central values of 
rric = 1-5 GeV and /i^ = rric, the hadronic contribution obtained with the GRV-HO 
set is 20% of the pointlike component, while the hadronic and pointlike contributions 
are comparable when using the LACl set. 

We also calculated charm electroproduction cross sections with the anti-tagging 
condition imposed, as in ref . [|10[ ; the results can be obtained from those presented in 
tables § and ^ by multiplying by a factor of 0.7, which is approximately independent 
of the parameters of the calculation. 

We have checked that varying the factorization scale of the photon in the pointlike 
component has a very small effect on the results. 
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mc = 1.2 GeV 


rric = 1.5 GeV 


= 1.8 GeV 




mc/2 


rric 


2mc 


mc/2 


rric 


2mc 


mc/2 


rric 


2mc 


GRV-HO 


0.878 


0.320 


0.176 


0.277 


0.129 


0.078 


0.115 


0.061 


0.039 


LACl 


6.697 


2.258 


1.204 


1.796 


0.774 


0.450 


0.646 


0.319 


0.196 



Table 6: Hadronic component of the total charm cross sections (^h) in ep 
collisions, at Ecm = 314 GeV, for various input parameters, as indicated. 
The MRSA set for the densities in the proton was used. 



In our analyses, we have kept the factorization scales iJ,p and /i^ fixed, since parton 
densities are usually available only for momentum scales above ~ 5 GeV^. There are 
actually two independent next-to-leading order parametrizations of densities in the 
proton valid for < m^: an extension of the MRSA set [^, which we will call 
MRSA modified, and the GRV-HO set . We have used these two parametrizations 
to perform a study of the factorization scale dependence of the charm cross section. 

The MRSA modified set is obtained by ordinary Altarelli-Parisi evolution of the 
MRSA set, supplemented with some specific assumptions on non-perturbative effects. 
On the other hand, the GRV-HO set has no non-perturbative assumptions. It is 
obtained by perturbative evolution, starting from "valence-like" inputs at Qq = 0.3 
GeV^. No attempt is made to mimic the \ow-Q^ data. The parametrization of non- 
perturbative effects inserted in the MRSA set is, strictly speaking, in contrast with 
the prescriptions of the factorization theorem. Nevertheless, we verified that these 
phenomenologically motivated inputs are numerically not important for factorization 
scales around the charm mass and larger. As far as the GRV set is concerned, the 
validity of the Altarelli-Parisi evolution equations and the perturbative expansion 
at the very low momentum scale at which the evolution is started is questionable. 
However, we believe that both sets are well suited for the purpose of estimating the 
uncertainties arising from the factorization scale variation. 

In fig. P we plot the ratio 

R = ^^^^ (4.3) 

= 2mJ 

as a function of ^. The dot-dashed curve is obtained with /x^ = fi, and keeping 
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Figure 9: The ratio R, defined in the text, for different choices of factor- 
ization and renormalization scales. MRSA (a) and GRV (h) densities in the 
proton are used. 



Hp = jjL^ = 2mc, the dashed curve has jdp = /j, /Jn = rric and /i-j, = 2mc, and the sohd 
curve has fip = fi^ = fi^ = /i. In the case of the MRSA parton densities, we also show 
(dotted curve) the same quantity, for iip = = = /i, with the parametrization 
of non-perturbative effects removed. In eq. ( [4.3|) a is the full (poinlike plus hadronic) 
cross section for an untagged electron and = 1.5 GeV. We use the GRV-HO set for 
the densities in the photon. When varying the hadronic factorization scale, the cross 
section, in the case of the MRSA modified set, decreases approaching the very low 
momentum region, where it increases again. This effect is essentially due to target- 
mass corrections in the parton densities, which sizeably harden the densities at low 
scales, and which are not present in the GRV set. The solid line is the result for the 
simultaneous variation of all the scales involved, including the photon factorization 
scale. We stress therefore that in this case we are probing also photon densities at 
low momenta. From this study we conclude that variation of the factorization scale 
brings about an additional uncertainty of 40% in the direction of lower cross sections. 

In tables ^ and ^ we present the pointlike and hadronic components of the total 
cross sections for bottom production at Eqm = 314 GeV. The defalut values corre- 
spond to all the scales equal to m;,, and = 4.75 GeV. As expected, in this case the 
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rub (GeV) 




DEF 


mb/2 


2mb 


mb/2 


2mb 


4.5 


5 


MRSA 


4.640 


5.332 


4.098 


4.391 


4.773 


5.541 


3.917 


CTEQ2MF 


4.814 


5.469 


4.283 


4.554 


4.958 


5.720 


4.079 


MRSD-' 


4.856 


5.580 


4.288 


4.649 


4.963 


5.839 


4.073 


CTEQ2ML 


5.123 


6.110 


4.418 


4.728 


5.290 


6.113 


4.325 



Table 7: Pointlike component of the total bottom cross sections (nb) in ep 
collisions, for various input parameters, as indicated, and -EcM = 314 Ge V. 













mb (GeV) 




DEF 


mb/2 


2mb 


mb/2 


2mb 


mb/2 


2mb 


4.5 


5 


GRV-HO 


1.115 


1.510 


0.851 


1.077 


1.146 


1.018 


1.177 


1.408 


0.892 


LACl 


2.842 


3.947 


2.142 


2.861 


2.840 


3.067 


2.675 


3.792 


2.157 



Table 8: Hadronic component of the total bottom cross sections (nb) in ep 
collisions, for various input parameters, as indicated, and Eqm = 314 GeV. 
The proton parton density set is MRSA. 



sensitivity to scale variations is much smaller than in the charm case. The sensitivity 
to the choice of parton densities is even less important than that arising from scale 
variation. This is because the small-x region, where the various sets mostly differ, is 
not deeply probed in the case of bottom production. The differences between results 
from various parton density sets come mainly from the different values of A. 

The hadronic component, displayed in table ||, is at most 75% of the pointlike 
component, even when the LACl set is used. Again, this is related to the fact that 
the small- a; region, where LACl gluon density is much larger than the GRV-HO one, 
is not so important in this case. In the case of bottom, electroproduction cross sections 
with the anti-tagging condition imposed on the scattered electron are approximately 
a factor 0.6 smaller than those presented in tables ^ and H, independently of all other 
parameters. 
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5. Conclusions 

Charm and bottom photoproduction cross sections are theoretically in a better 
shape than corresponding hadropro duct ion cross sections. We have presented results 
which indicate that the small-x effects should be under control at HERA energies, 
even for charm production, and that a measurement of the pointlike component of 
the charm photoproduction cross section at HERA could help to discriminate among 
various possible small-x behaviours of the gluon densities in the proton and in the 
photon. 

We found that the most important source of uncertainty is coming from the 
hadronic component of the cross section. If the gluon density in the photon is as 
large as the LACl parametrization suggests, then the hadronic component will be 
comparable or larger than the pointlike component for energies above 100 GeV for 
charm, and for energies of the order of 280 GeV for bottom. Hopefully, by looking 
at the structure of the events, experiments may be able to separate the pointlike and 
hadronic components, and therefore assess their relative magnitude. 

We also performed a study on total electroproduction cross sections in the Weiz- 
sacker- Williams approximation, including recent parton density parametrizations, the 
effect of an anti-tagging condition on the electron, and a more thorough illustration 
of the sensitivity to the various physical parameters. We find results that are in 
substantial agreement with earlier studies (ref. [TO). 
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